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ABSTRACT. The article provides an explicit algebraic expression for the generating 
function of walks on graphs. Its proof is based on the scattering theory for the differential 
Laplace operator on non-compact graphs. 



1. Introduction 

The concept of a generating function is known to be a very important tool in com- 
binatorics, probability, and number theory. Associated methods reduce the solution of 
combinatorial or probabilistic problems to the study of particular properties of the gener- 
ating function which can be performed by methods of function theory and analysis. For 
an introduction to this method the reader may consult the books [6], [8], l32l . A number 
of solved and still unsolved combinatorial problems, where the generating function plays 
a central role can be found in the article [25 ]. In the probabilistic context we mention the 
solution of the problem whether a simple random walk on Z d is recurrent or transitive 
by an analysis of the generating function (see, e.g., (9)). Some further examples will be 
discussed below in Sections|2]andEl 

The present work is devoted to the determination of the generating function for walks 
on graphs (both in combinatorial and probabilistic contexts). Walks on graphs are con- 
sidered, in particular, in GJ, 0, fT8ll . l26l Section 4.7], I33L In a custom setting random 
walks on graphs are defined as Markov chains on the vertices of the graph. The transi- 
tion probability from one vertex to another is assumed to be non-zero if and only if these 
vertices are adjacent. For a survey of the theory of random walks on graphs, see fTTll . 

We consider a slightly different but closely related model of random walks on graphs, 
where the states are chosen to be the edges of the graph. Transitions between different 
states are determined by stochastic (that is, Markov) matrices M(v ) prescribed at every 
vertex of the graph. The graphs are assumed to be non-compact, that is, besides a finite 
number of edges (or "internal lines" i S X) to have a non-empty set £ of "external lines" 
which serve as entries or exits for random walks. More precisely the model will be 
described in the following section. A relation between this model of random walks and 
random walks on vertices is explained in the Appendix below. 

Consider an arbitrary positive weight on the graph (that is, a map assigning to any 
edge i of the graph a positive number Oj). We will call a = {aj}j g j G a penalty 

vector. With j3 being a complex parameter we define a generating function T ee i ((3) of 
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walks from an external line e' 6 £ to an external line e S £ as 

= ^[MMW-^ [M^at-i)]^,^ . . . 

...[M{v{)\ i2M €-^ [M(«o)] ilie /, 

where the sum is taken over all walks {e , ii, «2j • • • >*JV-l)*JV) e } from e' to e, the set 
t>o, ... f at is the ordered list of vertices (with possible repetitions) visited by the walk, 
ik are the corresponding internal lines traversed during the walk (again with possible 
repetitions). In the context of the generating functions the weight exp{— fiai} can be 
viewed as a penalty factor for traversing the edge i during a walk. 

More generally, one can also consider penalty vectors depending on the direction in 
which a given edge is traversed by the walk. The corresponding generating function will 
be discussed in Section[6]below (see Theorem l6.7t . 

The main result of the present work (see Theorem 16.21 below) provides an explicit 
algebraic expression for the generating function of walks on graphs. Its proof is based 
on the scattering theory for the differential Laplace operator on non-compact graphs and 
the corresponding methods developed by the authors in lfT2l . lfl~3l . fl4l . (I5i . lfl6l . In the 
context of differential operators the weights a% will be interpreted as the metric lengths of 
the edges i. 

The generating function is determined by analytic continuation of the scattering matrix 
to complex values of the spectral parameter. This result is very reminiscent of a similar 
result in relativistic quantum field theory in the context of vacuum expectations of prod- 
ucts of quantum fields. The analytic continuation of the Wightman distributions 1011 to 
the Euclidean points (the so called Wick rotation 1 30 1) results in the Schwinger functions 
ll24l . Conversely, by a result [20|, [21 1 of K. Osterw alder and one of the authors (R. S.) the 
Schwinger functions give rise to Wightman distributions. In the bosonic case Symanzik 
and Nelson have shown that the Schwinger functions describe a stochastic theory (see 
0, 1191 . ll28ll . and references quoted there). 

We expect that the model of walks on graphs considered in the present article may be of 
interest in the context of optimization of traffic flows and in telecommunication networks, 
where the transition matrices M(v) determine a proportion of the traffic or signals to be 
transmitted in a given direction. 

The article is organized as follows. In Section |2] we will give definitions of walks on 
graphs and of associated generating functions. Also we present several examples relat- 
ing the generating function to combinatorics. In Section |3] we will revisit the scattering 
theory of differential Laplace operators on graphs. Section |4] is devoted to the proof of 
the combinatorial Fourier expansion formula (I4.18t . Theorem 14.21 proves the absolute 
convergence of the Fourier series and Theorem l4. lOl expresses the Fourier coefficients as 
sums over the walks on the graph. In Section|5]we will consider the analytic continuation 
of the scattering matrix with respect to the square root of the energy (that is, the spectral 
parameter). In Section[6]the generating function will be expressed in terms of the scatter- 
ing matrix for a Laplace operator with boundary conditions determined by the transition 
matrices M(v). In SectionQwe will turn to random walks on graphs. By means of the 
generating function we will calculate several mean values associated to this probabilistic 
set-up. 

There are several further models which can also be treated by the methods of the 
present work. In particular, choosing the matrices M(v) as independent random vari- 
ables one obtains a model of random walks in random environment. Further, the graph 
itself can be chosen to be random (see, e.g., [4]). We note that random graphs have been 
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used to model the spread of epidemics like AIDS, see, e.g., the article and further 
references quoted there. 



We consider a finite, connected and non-compact graph Q = (V,X, £, d), where V = 
V{Q) is a finite set of vertices, X is a finite set of internal lines, £ is a finite set of external 
lines. The elements of the set X U £ are called edges. The boundary operator d assigns to 
each internal line i £ Ian ordered pair (v\, v%) of vertices (possibly equal) and to each 
external line e G £ a single vertex v. The vertices v% := d~(i) and V2 ■= d + (i) are called 
the initial and terminal vertex of the internal line i, respectively. This obviously induces 
an orientation on each of the internal lines and this will become relevant below. 

The vertex v = d{e) is the initial vertex of the external line e. If d(i) = (v, v), then 
% is called a tadpole. To simplify the discussion, in what follows we will assume that the 
graph Q contains no tadpoles. 

Two vertices v and v' are called adjacent if there is an internal line i £ I such that 
either (v,v') = d(i) or (v',v) = d{i). A vertex v and the (internal or external) line 
j £ JU£ are incident if v G d(j). The degree deg(u) equals the number of (internal or 
external) lines incident with the vertex v. 

We do not require the map d : X — > V x V, £ — > V to be one-to-one. In particular, 
any two vertices are allowed to be adjacent to more than one internal line and two different 
external lines may be incident with the same vertex. 

Given an arbitrary vector a = {aj}j g x S M} 1 ^ with strictly positive components, we 
will endow the graph with the following metric structure. Any internal line i £l will be 
associated with an interval [0, a*] with ai > such that the initial vertex of i corresponds 
to x = and the terminal one - to x = ai. Any external line e € £ will be associated 
with a half-line [0, +oo). The number aj can be viewed as the length of the internal line 
i. 

A nontrivial walk w on the graph Q from e' G £ to e € £ is a sequence 



of edges such that 

(i) vq := d(e') G <9(ii), vat := d(e) G 9(ijv). and for any fc G {1, . . . , TV— 1} there 
is a vertex -u^ G V such that G <9(ifc) and G d(ik+i); 

(ii) u fc ^ v k+1 for all A; G {0, . . . , N - 1}. 

The number N is the combinatorial length |w| com b G N and the number 

N 



is the metric length of the walk w. 

Example 2.1. Let Q = (V,X,E,d) with V = {v ,vi}, X = {i}, £ = {e}, d(e) = v , 
and d(i) = (vq, vi). Then the sequence {e,i,e} is not a walk, whereas {e,i,i,e} is a walk 
from e to e. 

Proposition 2.2. Given an arbitrary nontrivial walk w = {e',ii, . . . e} there is a 
unique sequence {vk} k=0 of vertices such that vq = d(e') G d(i\), vn = d(e) G d(ijy), 
v k G d(ik), andvk G d(i k+ i). 



2. Walks on Graphs 



W,h, ■ ■ ■ ,iN,e} 
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Proof. Assume on the contrary that there are two different sequences {v k}k=o an( ^ Wk}k=\ 
satisfying the assumption of the proposition. This implies that there is a number K G 
{0, . . . , N - 2} such that v k = v' k for all k G {0, . . . , K} but v K +i / v' K+1 . Obvi- 
ously, the vertices vk, vk+i, and v' K+l are incident with the same edge. Thus, either 
vk = vk+i or vk = v 'k+1' which is a contradiction. □ 

We emphasize, that at any vertex of the sequence {v k } k=0 associated with a nontrivial 
walk w, the walk is either "reflected" or "transmitted". 

A trivial walk w on the graph Q from e' G £ to e G £ is the tuple {e',e} with 
d(e) = d(e'). Both the combinatorial and the metric length of a trivial walk are zero. 

A walk w = {e', i\, . . . , i^, e} traverses an internal line i G X if i k = i for some 
1 < k < N. It vw/to the vertex v if either w = d(e') or v = d(e) or v is incident with at 
least one internal line traversed by the walk w. 

The score n(w) of a walk w is the set {nj(w)}j g j with rij(w) > being the number 
of times the walk w traverses the internal line i G X. Any trivial walk has the score 
n = 0:={0,...,0}. 

Let W e , e ' = yVe,e'(G), e, e' G £ be the set of all walks wonfi from e' to e. In 
particular, the set W ei e' is infinite for all e, e' G £ if X / and the graph Q is connected. 
By reversing a walk w from e' to e into a walk w rev from e to e' we obtain a natural one- 
to-one correspondence between W e ^i and W e ' ie . Obviously, |w| = |w rcv | and n(w) = 
n(w rev ). 

Let S{v) C £ IJT denote the star graph of the vertex v G V, i.e., the set of the edges 
adjacent to v. Also, by S-(v) (respectively S + (v)) we denote the set of the edges for 
which v is the initial vertex (respectively terminal vertex). Obviously, S + (v)nS- (v) = 
since Q does not contain tadpoles by assumption. 

To every v G V we associate an arbitrary deg(v) x deg(v) matrix M(v) with com- 
plex entries [M(v)]j 1 j 2 , where ji,j2 G S(v) are edges incident with the vertex v. The 
collection of such matrices for all v G V will be denoted by M = {M{v)} v& v{g)- 

Now to each non-trivial walk w = {e', ii, . . . , in, e} from e' G £ to e G £ on the 
graph we associate a weight W(w) by 

|w| comb -l 

(2.1) W(w) = [M(^| Wcomb |)] e • J] [%)W.,H • [ M ^)]n,e' , 

I I comb , 

K=l 

where t> = <9(e'), ^| w | comb = ^(e), v k with G {1, . . . , |w| comb - 1} is the vertex 
incident with the internal line i k as well as the internal line i k +i- To a trivial walk w = 
{e',e} we associate the weight 

(2.2) W(w) = [M(3(e))] c>e , . 

Definition 2.3. The generating function of walks from e' G £ to e G £ on the graph Q 
associated with the collection M. = {M(v)} v< =v is defined as 

(2.3) Te.e'00) = Yl W(w)e-^ = ^(w)e- /3 ^ (w) ' a) , 

wew ee / wew ee / 

|w| = (n(w),a) := ^n,(w)aj. 

For given a walk w is called relevant if W(w) / 0. The set of relevant walks from 
e' to e is denoted by W e>e '(M). 
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Proposition 2.4. There is (3q > such that the series d2.3t converges for any e,e' £ £ 
and all (3 £ C w/?/j Re /3 > (3q. Moreover, 

(2.4) lim T<AP) = i mam '* " dle, = ale ' } ' 

Re/3too ' 10 otherwise. 

Definition 12. 31 suggests that we write W e>e ' as an infinite union of disjoint, non-empty 
sets by grouping together all walks w with the same score rt(w), 

W e , e '(s) = {w £ W e>e > | n(w) = n} 

such that 

(2-5) We,e' = (JWe,e>(n). 

n 

Note that these sets depend only on topology of the graph Q and are independent of its 
metric properties. Also if w £ W eie '(n) then w rcv £ W e ' )e (n)- W e , e '(Q) = if and 
only if 5(e) ^ d(e'). 

For the proof of Proposition |0] we need the following rather obvious fact: 

Lemma 2.5. The sets W eje '(n) are finite. Let 

fte fota/ number of internal lines traversed by any walk w £ W e>e /(n). 77je number 
of different walks in yV e ,e'(jk) satisfies the bound 



(2.6) |W e , e '(n)| < 



in ! 



n 



Set 

(2.7) r e , e ,(n)= £ W(w) 

if W e)e '(n) is nonempty and T eje /(n) = whenever W e ,e'(lL) = 0- Observe that T ejf >/(n) 
does not depend on the metric properties of the graph, i.e., is independent of the lengths 
of internal lines a. 

For given e, e' £ £ consider the set of scores of all walks from e! to e, 

(2.8) M e ,e' = {q I there is a walkw £ W e . e '(n)} . 

Since n(w) = n(w rcv ), we have N e ,e' = N e ',e- 

With this notation we have the following equivalent representation of (12.31) : 

(2.9) T ey (f3)= r e , e ,(n)e-^. 

Obviously, the series in (12.31) converges absolutely if and only if the series in (I2.3t does. 
Proof of Proposition \2.4\ Observe that 

|ni+i 



max ||M(u)|| ] 



|n|Re^a miE 



w£W„ e ,(n) weW„ e ,(n) 
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where 

(2.10) < a m i n := min a,. 

From Lemma l2~5l and using the identity 

Inl! 



(2.11) E TT~ = l X l Ar ' iVGN 

|n|=JV * eJ 



we, therefore, obtain 



E ( E W(w)e~^ 
neM , weW e ./(n) 



< E (max||M(V)|| J e -^/3a min | T |JV_ 



JV=0 

This series converges for all [3 E C with Re /3 > /3q, where 
(2.12) fa > — ( maxlog ||M(u)|| + log |X 

We mention also the following simple result: 



□ 



Lemma 2.6 (Time Reversal Invariance). If all matrices M(v) are symmetric then so is 
the matrix T(f3) with matrix elements T ee i(f3) for all large Re/5 > 0. If all M(v) are 
self-adjoint, then so is T{(5)for all large (3 > 0. 

Definition 2.7. The family of matrices A4 is called combinatorial if every matrix entry of 
every matrix M(y) equals either zero or one. 

If M. is combinatorial, the weight W(w) of an arbitrary walk w is either zero or one 
and we have the following simple result. 

Lemma 2.8. If M is combinatorial and W eje ' (M) finite then 

Te,A0) = \We,AM)\, 
i.e., the number of relevant walks from e' £ £ to e E 8. 

We now provide some examples, which relate our formulation to well known com- 
binatorial contexts. Viewing 1? as a subset of K 2 , for an arbitrary n E N consider the 
set 

V n = {(x u x 2 ) E Z 2 \ < x 2 < xi < n) . 

We consider the non-compact graph Q n = (V n ,I,£,d), where £ = {e, e'}, d(e') = 
(0, 0) and d(e) = (n, n), and the vertices v\ E V n and v% E V n are adjacent if and only 
if the Euclidean distance between these vertices is not larger than y2, I ^i — v 2 \ < \/2. 
Therefore, the set of internal lines X consists of of all intervals joining the points of V n and 
having Euclidean distance not greater than y/2 (see Fig.Q. The metric distance between 
two adjacent vertices will be assumed to be equal 1, that is, a% = 1 for all i EX. 

Example 2.9 (The Catalan numbers). The number 

1 (in 



C i — r „ I , "ell 
n + 1 V n 
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is called the (n — l)-th Catalan number (see, e.g., [29| and pp. 219 — 229 in \21\). Set 
/C CataJan = {(1, 0), (0, 1)}. For an arbitrary vertex v G V n of the graph Q n and arbitrary 
j G £ fl 2 adjacent to the vertex v we set 

v' — v G Z 2 if j £ I and is adjacent to the vertex v', 
Xv{j) = Ul,0)eZ 2 ifj = e, 
,K0)£Z 2 ifj = e'. 

Let K, Catahn = {(1, 0), (0, 1)} C 1? and 

r M c^n ( v)] _ f 1 ifXv(h) e K Catalan and - X v{h) G 
L 1 | otherwise. 

The setyV e: e'{-M. Catalan ) is, obviously, finite. Therefore, the generating function T ee i (/?) 
jj entire. For given n G N number T e e j (0) w f/je (n — Catalan number. The 
three other matrix elements T e ^((3), T e ^ e (/3), and T e / ie /(/3) vanish identically. 

In the next example we continue with the same notation. 

Example 2.10 (The Schroder numbers). The Schroder numbers ( see, e.g., [ 29 1 and p. 1 78 
in 1271 ) can be defined by the recurrence relation 

n-1 

S n = S n -i + SkS n -k-i with So = 1. 

fc=0 

Le? /C Sd,r6der = {(1, 0), (0, 1), (1, 1)} D K, Catalan and 

r Schroder 



\M Scbr6der (v)] . . 



I 1 ifXvtii) G /C™- anrf - Xu (j 2 ) G /C™ ej ', 
JU2 I o otherwise. 



Obviously, W e)e >(M Catalan ) C W eie /(A^ Sciir ° der ) is again a finite set. For given n G N 
number T e e < (0) is now n-f/j Schroder number. The three other matrix elements 
T e<e (f3), T e ^ e (P), and T e ' j£ i ((3) vanish identically. 
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For the next two examples consider the sets 

V+ = {(xi,x 2 ) G I?\ < xi < n, < x 2 < n) , n G N. 

Let = (V+,1, £,d) be the non-compact graph with £ = {e, e'}, d(e') = (0,0) 
and <9(e) = (n, 0), and the vertices v\ G V and v 2 G V are adjacent if and only if 
the Euclidean distance between these vertices is not larger than ^/2, \v\ — v 2 \ < \/2. 
Therefore, the set of internal lines X consists of of all intervals joining the points of 
and having Euclidean length not greater than \[2 (see Fig. EJ- 

Example 2.11 (Dyck paths). Let K Dyck = {(1, 1), (1, -1)} and 

lM Dyc k(v)] _ fl tfXvUl) e K2 Dyck and - Xv (j 2 ) G JC Dyck , 
L %,i 2 | otherwise 

if neither j% nor j 2 are external lines. We set 

'l ifXv(ji)elC D y ck , 



i MDyCk (v)] n , e , 



otherwise 



if j 2 = e' and 



i if - xvih) g ic Dyck , 



\M Dyck {v)} . -- , 

ej2 otherwise 



if h = e. 

Obviously, W e ^ e '(M Dyck ) is a finite set. Therefore, T efi i{(3) is entire, T eje /(0) is the 
number of Dyck paths on the graph Q+. A discussion of Dyck paths can be found in\\ \ \. 

Example 2.12 (Motzkin numbers). The non-compact graphs are the same as for Dyck 
paths in Example 12. 1 11 The Motzkin numbers (see, e.g., Q and Problem 6.37 in [27]) 
can be defined by the recurrence relation 

n-2 

M n = M n _i + M fc M n _ fc _ 2 with M = M l = 1. 

k=0 
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Set K, Motzkin = {(1, 1), (1, -1), (1, 0)} D KPy ck and 

\M Mo '^(v)] = { 1 ifXv{3l) G KM ° tZkin aM ~ Xv{n) ^ /CM ° tZfan ' 
L >ih,h | otherwise. 

Again the set W e!e '(Ai Mo ") is finite and, therefore, T e ^i{(3) is entire. For given 
n E N f/ie number T ee i (0) is f/je n-f/j Motzkin number. 

3. Laplace Operators on Graphs 

In this section we will recall the theory of Laplace operators on a metric graph Q and 
the resulting scattering theory (see [12), llT3l . fl4l . lfT5l . (T6'| for further details). 

Given a finite non-compact graph Q = (V,T, £, d) with a metric structure a = {aj}j g j 
consider the Hilbert space 

(3.1) n = n(£,i,a) = n £ en I , n £ = ($n e , h x = Qtu, 

ee£ iel 

where H e = L 2 ([0, oo)) for all e E £ and Hi = L 2 ([0, a*]) for all % E J. By £>, with 
j E £ U X denote the set of all tpj E Hj such that ipj(x) and its derivative ipj(x) are 
absolutely continuous and ip"(x) is square integrable. Let denote the set of those 
elements tpj of T>j which satisfy 

^■(0)=^(a j )=^(0)=^(a i )=0 for j el 

and 

^■(0) = ^(0) = for j€f. 
Let A be the differential operator 

(3.2) (A» . (x) = j ETU £ 
with -0 = {V'ilieJu^ in the domain 

D° = V] c W. 

It is straightforward to verify that A is a closed symmetric operator with deficiency 
indices equal to \S\ + 2|X|. 

We introduce an auxiliary finite-dimensional Hilbert space 

(3.3) K = K{£a)=K £ ®1C { I ) 

with ]C £ = C' £ and /Cj^ = C' 1 '. The subspaces K,j ^ we associate with initial vertices 
of the internal lines i EX, the subspaces ACj + ^ with the terminal vertices. Let d K denote 
the "double" of K, that is, d K = K@lC. 
For any ip E V := Vj we set 

(3.4) :=±®f £ d JC, 
with 

/{^e(0)} ee A / {^(0)} e6 £ \ 

(3-5) ±= M(0)W f= {^(0)}*ez G/C. 

\{^i(ai)}iex/ \{-^( a i)}ie2:/ 
Here the vector notation is used with respect to the orthogonal decomposition (I3.3I) . 



to 
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To define the Laplace operator on the graph Q consider the family tp = {ipj}j£Eui of 
complex valued functions defined on [0, oo) if j £ £ and on [0, a.j] if j £ X. Formally the 
(self-adjoint) Laplace operator is defined as 

d 2 

(3.6) (A(A,B,a)?P) j (x) = ^iP 3 (x), jelUS 

with the boundary conditions 

(3.7) Axjj_ + B^_ = 0. 

By definition A and B are any complex (|£| + 2\I\) x (|£| + 2\I\) matrices such that 

(i) the matrix (A, B) has maximal rank, 

(3.8) , 

(ii) the matrix AB' is self-adjoint. 

Here and in what follows (A,B) will denote the + 2\1\) x 2(\£\ + 2\1\) matrix, 
where A and B are put next to each other. 

The scattering matrix S(k) = S(k; A, B,a) associated to A(A, B,a) has the follow- 
ing interpretation in terms of the solutions to the Schrodinger equation (see [ 12 1 and 
[15 1). Consider the solutions ip k (k) (k £ £) of the stationary Schrodinger equation for 
-A(A, B,a) at energy k 2 > 0, 

-A(A,B,a)ij k (k) = kV(k) 



of the form 



S(k) jk e lkx for je£,j^k 

e- lkx + S{k) kk e lkx for j e£,j = k 



(3.9) rf(x;k) 

' a(k) jk e lkx + (3(k)jke~ lkx for j £ X. 

Thus, the number S(k)j k for j ^ k is the transmission amplitude from channel k £ £ to 
channel j £ £ and S^k)^ is the reflection amplitude in channel k £ £. Their absolute 
squares may be interpreted as transmission and reflection probabilities, respectively. The 
"interior" amplitudes 

a(k) jk = a(k;A,B,a) jk , I3(k) jk = (3(k;A,B,a) jk 

are also of interest, since they describe how an incoming wave moves through a graph 
before it is scattered into an outgoing channel. 

The condition for the ijj k {E) (k £ £ ) to satisfy the boundary conditions i3.lt imme- 
diately leads to the following solution for the scattering matrix S(k) : K,£ — ► K.£ and 
the operators a(k) and /3(k) acting from Ks to Kj. Indeed, by combining these operators 

into a map Ks to K, = Kg © K,^ © we obtain the linear equation 

fS(k)\ fP 

(3.10) Z(k;A,B,a) a(k) = -(A-ikB) 

V/3( k )/ \°y 

with 

(3.11) Z(k;A,B,a) = AX (k; a) + ikBY (k; a) 
where 



(3.12) X{k;a) = I I , F(k;a) 










I 


—I 


e ika 


e -ika 
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±ika 



The diagonal |X| x |X| matrices e - are given by 
(3.13) ej lka = <5 jfc e ±lk ^ for j, k € I. 

Theorem 3.1 (= Theorem 3.2 in O). For any k G M 

Ran(A-ikB) I I C Ran Z(k; A, B, a). 



77ims, equation (I3.10t a solution even if (let Z(k; A, B, a) = Ofor some k £ 1. 77z/s 
solution defines the scattering matrix uniquely. Moreover, 

T 

(3.14) S(k) = -(I 0)Z(k;4,B,a)- 1 P K i er , (MAfl) ( J 4-ikS) [o 



,0, 



« unitary for all k 6R \ {0}. 



In the case with no internal lines (X = 0) the relation (I3.14t for the scattering matrix 
simplifies to 

(3.15) S(k; A, B) = - (A + ikB)- 1 (A - ikB) . 

Proposition 3.2. 7/'det(A + ikB) = Ofor some k G C, then k = ik with k G M. For any 
sufficiently large p > ?/zere w a constant C p > smc/j ?/za? 

(3-16) IKA + ikS)- 1 !! < ^(l + iki)- 1 

/or a// k G C w/f/j | k| > p. 

Proof. Assume that det(A + ikB) = for some k G C with Re k ^ 0. Then also 

det(A t - ikB*) = det(A + ikB) = 0. 
Therefore, there is a \ 7^ such that 
(3.17) (A* - ikB+)x = 0. 

In particular, we have (BAt - ikBBt) x = 0. Therefore, since BA^ is self-adjoint, we 
get 

(X,BA^ X ) = (X,BB^ X ) Imk, 
( X ,BBt x ) Rek = 0. 

The second equality implies that x G KerB^. Then, by (I3.17I) . % £ Ker A'. Since the 
matrix (A, B) is of maximal rank, we have Ker A^ n Ker B^ = {0}. Thus, x = which 
contradicts the assumption and, hence, Re k = 0. 

Since det(A + ikB) is a polynomial in k, it has a finite number of zeroes. Take an 
arbitrary p > such that all its zeroes lie in the disk |k| < p. Using the matrix inverse 
formula we represent any element of (A + ikB) -1 as a quotient of two polynomials of 
degrees \£\ + 2|X| — 1 and \£\ + 2|X|, respectively. In turn, this implies the estimate 
(13161 . □ 

Theorem 3.3. The scattering matrix S(k) = S(k; A, B,a) is a meromorphic function in 
the complex k-plane. In upper half-plane Im k > it has at most a finite number of poles 
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which are located on the imaginary semiaxis Re k = 0. Outside these poles the scattering 
matrix is holomorphic for all Im k > and determined by the relation 

(3.18) 5(k)=-(l 0) Z{k-A,B,a}~ 1 {A-ikB) 

Proof. Assume that det Z(k; A, B,a) = for some k G C with Im k > and Re k / 0. 
This implies that the homogeneous equation 



Z(k;A,B,a) [a = 




has a nontrivial solution with s G Ks and a, f3 G C'^L Consider the function ip(x) 
{ipj(x)}jelu£ defined by 



lj)j(x) 



Sje lkx for j G £, 

aje lkx + f3 je - lkx for j G 2. 



Obviously, ip(x) satisfies the boundary conditions J3.7t . Moreover, ip G £ 2 (<?) since 
Im k > 0. Hence, k 2 G C with Im k 2 ^ is an eigenvalue of the operator A(^4, -B, a) 
which contradicts the self-adjointness of A(^4, B,a). 

Since det Z(k; A, B, a) is an entire function in k which does not vanish identically, 
from (I3.10t it follows that the scattering matrix 5(k) is a meromorphic function in the 
complex k-plane. To prove that the scattering matrix S(k) has at most a finite number 
of poles on the imaginary semiaxis {k G C| Re k = 0, Im k > 0} it suffices to show 
that the determinant det Z(k; A, B, a) does not vanish for all sufficiently large Im k > 0. 
To see this we set k = ik with k > and assume there is an unbounded non-decreasing 
sequence {K/JfceN such that 

detZ(iK k ;A,B,a) = for all k G N. 

Therefore, there is a sequence {xfcjfceN of normalized elements \k G /C such that 

X(iK k ;a)^A^Xk = K k Y(iK k ; a) ] B ] Xk- 
It is straightforward to verify that X(in k ; a) is invertible and 





-l 



R k := (^X (iK k ; a,y^j Y(iK k ;a)^ = | coth(Ka) — [sinh(Ka)] 

\0 — [sinh(/ta)] _1 coth(Ka) 

with a notation analogous to (I3.13I) . Thus, 

(3.19) - K k B^) Xk = n k {R k - l)B^ Xk 

for all k G N. Observe that \\R k — I|| = 0(e~ CKk ) for some c > as k — > oo. By 
Proposition l3.2l the operator A< — kB^ is invertible for all sufficiently large k. Moreover, 
\\(A^ - kB^^W < C with C > for all sufficiently large k. Thus, equation d3~T9l 
implies that \k — * which contradicts the assumption ||xfe|| = 1- ^ 

In the lower half -plane Im k < the scattering matrix may have poles with Re k ^ 
(see, e.g., Example 3.2 in El )- These poles correspond to resonances. 

The notion of local boundary conditions has been introduced in our article [12| and 
is discussed in more details in [15] and [16|. Local boundary conditions couple only 
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those boundary values of tp and of its derivative ip' which belong to the same vertex. The 
precise definition is as follows. 

With respect to the orthogonal decomposition K, = Kg © ICj ^ © K,^ any element z 
of /C can be represented as a vector 

/ {Ze}ee£ \ 

(3.20) z= • 

Consider the orthogonal decomposition 

/C = C v 

with C v being the linear subspace of dimension deg(v) spanned by those elements d3.20t 
of K, which satisfy 

z e = if e £ £ is not incident with the vertex v, 

(3.21) z\ = if v is not an initial vertex of 

z^ =0 if v is not a terminal vertex of i G X. 
Set d C v :=C V ®C V ^ C 2Ac ^ v ). By the First Theorem of Graph Theory we have 

deg(» = \£\ + 2\T\ 

v€V(G) 

such that 

© d C v = d K. 

vev(g) 

Definition 3.4. Given the graph Q = Q(V,2,£,d), the boundary conditions (A,B) 
satisfying (I3.8t are called local on Q if and only if there is an invertible map C : K, — > K, 
and linear transformations A(v) and B(v) in C v such that the direct sum decompositions 

(3.22) CA = Q)A{v) and CB = Q)B{v) 

hold simultaneously. Otherwise the boundary conditions are called non-local. 

For instance, for a single- vertex graph any boundary conditions are local. The bound- 
ary conditions considered in Example 3.4 of [ 16] are non-local. 

4. Combinatorial Fourier Expansion of the Scattering Matrix 

In this section we will perform a harmonic analysis of the scattering matrix with re- 
spect to the lengths a = {ajjiex £ of the internal lines of the graph Q. The 
main results of this section are presented in Theorems 14. 21 and 14. 1 01 In Theorem 14.21 the 
absolute convergence of the Fourier series for the scattering matrix is proved. Theorem 
14. 101 expresses its Fourier coefficients as sums over the walks on the graph. Combining 
these two results proves the combinatorial Fourier expansion formula (I4.18t . 

Throughout the whole section we will assume that the (topological) graph Q as well as 
the boundary conditions (A, B) are fixed. To carry out the analysis we will now treat a 
as a parameter which may belong to l^l or even C^L 

We start with the following simple but important observation. 
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Lemma 4.1. For arbitrary k > the scattering matrix S(k; A, B,a) is uniquely defined 
as a solution of (13.101) for all a G R™. Moreover, the scattering matrix is periodic with 
respect to a, 

S ^k; A,B,a + ^l) = S(k\ A, B,a) 

for arbitrary £ E Z^L 

Proof. It suffices to consider those a G for which det Z(k; A, B,a) = 0, since the 
claim is obvious when the determinant is non- vanishing. For a G (M_|_)' :r ' the fact that 
S(k; A, B,a) is uniquely defined as a solution of (I3.10t is guaranteed by Theorem 13. II 
The case of arbitrary a G R' 1 ' can be treated exactly in the same way (see the proof of 
Theorem 3.2 in fT2ln . 

The periodicity follows immediately from (13.101) and the fact that the matrices X(k; a) 
and Y(k; a) in dTTH are -periodic. □ 

Lemma 14.11 suggests to consider a Fourier expansion of the scattering matrix. The 
following theorem ensures the absolute convergence of the corresponding Fourier series. 

Theorem 4.2. Let k > be arbitrary. For all a G R'" 1 ' the Fourier expansion of the 
scattering matrix 

(4.1) S(k;A,B,a)= £ 4(k;i,B)e lkM 

nezm 

with 

(4.2) Sn(k;A,B) = Q^-Y J da S (k; A, B , a) e~ lk ^ 

[0,27T/k] 1^1 

converges absolutely and uniformly on compact subsets of R^L The Fourier coefficients 
( 14.21) vanish for all n = {rii}i e j G Z^l w/?/z nj < 0/or af feasf one i El. 

For the proof we need a couple of auxiliary results. Set 

A = {a = {ai} i( zj \ Reai G R, Ima.; > 0} C C' J L 

Lemma 4.3. For a/ry k > ?/ze determinant det Z(k; A, B,a) has no zeroes for all 
a G A. 

Proof. Assume there is a G A such that det Z(k; A, B,a) = 0. Then there are s G 
and a, (3 G C' 1 ' such that 

Z(k:A. B.a) I a I =0. 



Equivalently this gives 



(,1 + lki?) ( n I + (.1 - ikB) I [3 \ U. 

~ lka -f3 J \e lk ^a, 

The operator (A + ikB)^ 1 (A — ikB) is unitary for all k > (see the proof of Theorem 2.1 
in lU2l ). Since unitary transformations preserve the canonical Hilbert norm on C^"*" 2 ! 1 !, 
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we have 

INI 2 + N^ 1 " e" 2kImai ) + IA| 2 (e 2kIma! - 1) = 0, 

which implies s = and a = (3 = 0. □ 

Proposition 4.4. Let k > be arbitrary. For all 

a G clos(.A) := {a = {aj}j G j| Rea^ G R, Imaj > 0} 

scattering matrix S(k; A,B,a) is uniquely defined as a solution of (I3.10t and satisfies 
the bound 

(4.3) \\S(k;A,B,a)\\ < 1. 

Moreover, it is a rational function oft = {tjjjgi w/?/j ij := e lka % i.e. a quotient ofB{K£)- 
valued polynomials in the variables t{. Thus, for all a G clos(^l) the scattering matrix is 
-periodic, 

S (k;A,B,a + Y&) = S(k,A,B,a), leZ^. 

Proof. By Lemma FOl equation (I3.10t has a unique solution for all a G A. Equations 
(I3.11t and d3.12t imply that Z(k; A, B,q) is a polynomial function of the components 
of t. Obviously, Z (k; A, B , a)^ 1 is also a rational function of t. Thus, by (I3.10t the 
scattering matrix S(k; A, B,a) is a rational function of t. Thus, it is ^Z^-periodic. 
Using (13. 10b it is easy to check that this solution satisfies the relation 

S(k;A,B,a) \ I I 

a(k;A,B,a) \ = - (A + ikB)- 1 (A - ikB) \ f3(k;A,B,a) 
K e~ lk ^f3(k; A,B,a)J \e lk ^a(k; A, B, a) 

Since (A + \kB)^ 1 {A — lkB) is unitary we obtain 

S{k; A, B,a)^S(k; A, B,a) + a(k; A, B,a)^(I - e - 2kIm ^)a(k; A, B,a) 

i4A) +P(k;A,B,ay(e 2klm ^-I)P(k;A,B,a)=I, 

where Ima = {Imaj}j e x- From (14.41) it follows immediately that 

< S(k; A, B,a) J( S(k; A, B,a) < I 

in the operator sense. This proves the bound (I4.3t for all a G A. Recalling Lemma |4~T1 
completes the proof. □ 

A priori it is not clear whether the boundary values of the scattering matrix 5(k; A, B, a) 
with a G A coincide with those given by equation d3.14t for all a G W 1 '. The following 
lemma shows the "non-tangential continuity" of the scattering matrix S(k; A, B,a) with 
respect to a G clos(.A). 

Lemma 4.5. Let a G R' 2 ' and k > be arbitrary. For any sequence {aj}j^, a,j G A 
converging to a G R' 2 ' the relation 

(4.5) lim S(k; A, B, a,) = S(k; A,B,a) 

holds. 

For the proof we need the following elementary result. 
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Lemma 4.6. Let T n be a sequence of invertible operators on the finite-dimensional 
Hilbert space converging to the operator T. Then 

i^o^Ker T T n ^KerTt = P KerT T ^KerTt' 

where Pg denotes the orthogonal projection onto orthogonal complement in f) of the 
subspace £ C fj. 

Proof. Consider the operators T n and T as maps from (KerT) ± to (KerT^) ± . Since 
these maps are invertible, the claim follows from the obvious relation 

T~ l = T^ [i + T-^Tn-T)}- 1 . 



□ 



Proof of Lemma W3\ Introduce the shorthand notation 

Z(a) = Z{k;A,B,a) and S(a) = S(k; A, B, a) 
From Theorem B.ll and Lemma l4~3l it follows that 

'r 

S(a) = — (I 0) Z{a)- l P^ cvZ{ ^{A - xkB) \ 



,0, 



and 



S(a j ) = - (I 0) Z{a j )~ 1 {A - ikB) |o 

= -(1 0)Z(a j r l P^ erZ(&) ,(A-ikB)\0\. 
Thus, to prove the claim it suffices to show that 

(4.6) Urn (I 0)Z(a j r 1 P^ crZ{ t = (l 0) Z^P^ z( t . 

From Theorem 3.1 in (12 1 it follows that all elements z of Ker Z(a) satisfy Vgz = 0, 
where Vs is the orthogonal projection in K, onto Ks- Thus, 

(4.7) (I 0) Z^.)- 1 ^^ = (I 0)Pi eiZ{ a ) Z(a j r 1 Pi eTZ(A)t 
for any j G N and 

(4.8) (I 0) ZCa)- 1 ^^ = (I 0)^ z(a) Z(a)- 1 i^ rZ(a)t . 
By Lemma l4~6l we have 

^^KerZ(a) 2 fei) 1 P Ker Z (a)t = P Kev Z (a) Z (^) lp Ker Z(a)f- 
3 >0 ° 

Combining this with (14.71) and d4.8t we obtain (14.6b . □ 
For fixed k > 0, consider 

(4.9) F(t):=S(k;i,B,a) with £ = e^ k . 



Recall that by (I3.10t - (13.121) the scattering matrix S(k;A, B,a) depends on a only 
through t = e 1 - . The map a i— ► e'- k maps the set 

| a G C' 1 ' | a = {aj}j e i with < Re a, < 27r/k and Im > for all i G Z j . 
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bijectively onto the polydisc H)^ = {( G C| |£| < 1} ' :r ' ■ The interval (0, 2ir/k] is mapped 
onto the torus Tl x l = {C G C| |C| = 

Lemma 4.7. The function F belongs to the Hardy class H p (T^)for all p G (0, oo] gotc? 

Remark 4.8. We recall that an operator-valued function on a polydisc U> d is said to be 
inner if it is holomorphic in B d and takes unitary values for almost all points of T d C 
<9(B d ) (the so called distinguished boundary ofH d 1101 ). For d = 1 matrix-valued inner 
functions are studied, e.g., in [22]. In particular, an analog of the canonical factorization 
theorem for matrix-valued inner functions has been proven there. 

Proof. From Proposition l4.4l it follows that F is holomorphic in the punctured open poly- 
disc B^l \ {0}. By (03) we have \\F(t)\\ < 1 for all t G B |J| \ {0}. Therefore, the 
Laurent expansion of F contains no terms with negative powers. Thus, F is holomorphic 
inB^I. 

The bound (14.31) also implies that 

sup / \\F(rt)\\ P dm < 

re[0,l) Jjm 

for any p G (0, oo), where [i stands for the Haar measure on the torus T^l and 

sup sup ||.F(ri)|| < 1. 

For every t G T^l the operator F(t) is unitary, which means that F is an inner function. 

□ 

Proof of Theorem W2\ Since by Proposition l4.4I Fft) is a rational B(K,g) -valued function, 

it can be analytically continued as a meromorphic function on all of t G C' 1 '. Moreover, 

\x\ 

it is holomorphic in the polydisc B^ = {( G C| \(\ < 1 + s} for some e > 0. To show 
this, by Hartogs' theorem it suffices to consider the analytic continuation with respect to 
a single variable tj G C keeping all other variables fixed. By the bound (I4.3t all possible 
poles of this continuation lie outside a disc {tj G C| < r} with r > 1. 

In turn, this implies (see, e.g., Theorem 2.4.5 in [ 10 1) that the Taylor series of the 
function F(t) converges absolutely and uniformly for all t G T^L Combining this with 
Lemma POl proves the absolute and uniform convergence of the Fourier expansion d4.lt . 

By Lemma |4~71 the Fourier coefficients d4.2t satisfy S n (k; A, B) = for any n G Zl x l 
with rii < for at least one i G X. □ 

Definition 4.9. Given a non-compact graph Q = (V,I,£,d) to any vertex v G V = 
V(Q) we associate the single-vertex graph Q v = ({v},l v ,£ v ,d v ) with the following 
properties 

(i) Fv = 0, 

(ii) d v (e) = vfor all e G £ v , 

(hi) \£ v \ = degg(f ), the degree of the vertex v in the graph Q, 

(iv) there is an injective map fy v : £ v — > £L)1 such that v G do^ v (e)for all e G £ v . 

Since the boundary conditions are assumed to be local (see Defmition l3.4l) . we can con- 
sider the Laplace operator A(A V , B v ) on L 2 (Q V ) associated with the boundary conditions 
(A v , B v ) induced by (^4, B), see d3.22l) . By d3.15t the scattering matrix for A(A V , B v ) is 
given by 

S v (k) = -(A v + ikB v )- 1 (A v -ikB v ). 
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Now to each walk w = {e', i\, ■ ■ ■ ,in, e) from e! G £ to e G £ on the graph similar 
to (I2.lt we associate a weight VF(w; k) by 

(4. 10) W(w; k) = e lk ^ (w) ' a) VF(w; k) 
with 

l w lcomb 

(4.11) W(w;k) = [J 5„ fc (k) e(+)e( _). 

fc=0 

Here ei G are defined as 

(-) = J*« fc 1 (*ft)> if 1 < ^ < |w|comb, 

and 

(+) _ l^Vkdk+l)' if < ^ < l w !comb - 1, 

6fc [^(e'), if fc = |w| comb + 1, 

where the map ty v is defined by Definition 14.91 Note that W(w; k) is independent of 
the metric properties of the graph. Obviously, for a trivial walk w = {e',e} we have 
W(w; k) = ^(k) f -i (e)itt -i (e))) where v = 8(e) = 8{e'). 

Theorem 4.10. The matrix elements of the n-th Fourier coefficients (14. 2t are given by the 
sum over the walks with score n, 

(4.12) [S Jk (k;A,B)] eie >= Yl W(vs:k) 

weW ee /(n) 

ifW e ,e'(lk) is nonempty and [5„(k; A, -B)] e . e ' = whenever W e ,e'(ll) = 0- 

Proof. Obviously, it suffices to show that the n-th coefficient of the multi-dimensional 
Taylor expansion of the scattering matrix S(k; A, B, a) with respect to t = {ti}i e x G Bl x l 
with ti := e lkcii coincides with the r.h.s. of (14.1 21) . Recall that by Theorem l3.1l and Lemma 
I4.3l for all a G A the scattering matrix is given by 



(4.13) S(k; A, B,a) = — (l 0) {AX(k; a) + ikBY(k; a))" 1 {A - lkB) I I , 



,0, 



where X(k; a) and Y(k;a) were defined in (I3.12t . Obviously, 
(4.14) 

AX(k;a) + ikBY(k;a) = (A + ikB)£/(k; a) + (A - ikB)R(k; a) 

= {A + ikB) [I + (A + ikB) _1 (4 - ikB)12(k;a)?7(k;o)- 1 ]C/'(k;a), 

where 



(4.15) U(k;a) := [ I 

and 



e" lka 



'0 0^ 

R(k;a) :=X(k;a)-U(k;a)= | I 

,0 e lk ^ 



RANDOM WALKS ON GRAPHS 



19 



with respect to the orthogonal decomposition (13-31) - Equation (14.141) implies that 

(AX{k;a) + ikBY(k;a))~ 1 

oo 

= Ui^ay 1 ^ [-{A + ikB)- 1 {A - ikB)GH{W,a)] n {A + ikB)- 1 

?1=0 

with 

(4.16) C; - I I J and II (k: a) 










,ika 








gika 



such that R(k; a)f/(k; a) 1 = G-ff(k;a). Combining this representation with (I4.13t we 
obtain 

S(k;A,B,a) = 

(4 17) 00 

= J^(I 0)[S{k;A,B)GH(k;a)} n S{k;A,B) \0 

n=0 \ . 



where S(k; A, B) is defined by (I3.15I) . By the unitarity of S(k; A, B), the series converges 
absolutely for all a G A. 
Recall that 

S{k;A,B) = S{k;CA,CB) 
for every invertible C. It follows directly from Definition l3.4l that 

S(k;A,B)= S(k;A(«),B(t;)). 
veV{G) 



Plugging this equality in (I4.17t proves the claim. □ 

Remark 4.11. Theorem 14. 1 01 implies that the scattering matrix of the graph Q is deter- 
mined by the scattering matrices associated with all its single vertex subgraphs. This 
result can also be obtained by applying the factorization formula [ 14 1. 

Combining Theorems I4.2l and l4.10l we immediately obtain 

Corollary 4.12. Let a € (M-j-)'- 1 ' be arbitrary. For all k > the scattering matrix 
S(k; A,B,a) associated with the Laplacian A(A, B,a) on the graph Q has an absolutely 
convergent expansion in the form 

(4.18) S(k;A,B,a) e;el = £ W(w; k) = £ W(w;k)e lk H. 

Since W(w; k) is independent of the metric properties of the graph, it is natural to 
ev 

r e 

IP 



view (14.181) as a combinatorial Fourier expansion of the scattering matrix S(k; A, B,a). 
We will show that (I4.18t actually coincides with the Fourier expansion (14.11) in Theorem 



5. Analytic Continuation of the Scattering Matrix 

Recall that the scattering matrix S(k; A, B,a) is analytic in k for all Re k > and 
Im k > 0. In this section we will show that representation (14.181) for the scattering matrix 
can be extended to the complex plane. 
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Lemma 5.1. There is f3g > such that the series 

(5.1) [Sn(k;A,B)] e , e ,e lk ^ 

converges absolutely for all k G C wiY/z Re k > and Im k > /3o- Therefore, 
S(k;A,B,a) e , e ,= ]T [5„(k; A, 5)] e , e ,e lk ^ 

is a holomorphic function for all such k G C. 

Proof. From Proposition 13.21 it follows that for all sufficiently large j3 > there is a 
constant > such that the estimate 

holds for all k G C with Im k > (3, any v £ V, and any ei, e2 G ^ (see Definition I4.9I) . 
Therefore, for an arbitrary walk w G W eje ' (n) we obtain 

rM + l 



|W(w;k)| <C 







Thus, from (I4.12t it follows that 



\[Sn(k;A,B)} ete/ \<C^ +1 \W e , e/ (n)\. 
Therefore, from Lemma 1231 using the identity i2. lit we obtain the estimate 

(5-2) l[4(M,B)]e,e'l <C$ +1 \1\ 

\n\=N 

Recalling the definition ( I2.10t for a m j n estimate ( I5.2t implies that the series 

£ |[4(k;^l,i?)] eie ,e lk ^| < J] |[4(k;^,B)] eie ,|e-^ Imka ^ 

oo 

n)e,e'\ 



< ^ e -ATImka min ^ | 5 ( k; 

\n\=N 



N=0 R^K,e' 



converges for all k G C with 



bnk>#L loglC^}. 

This proves the claim with (3q = max{/3, □ 

The following statement is the main result of this section. 
Theorem 5.2. There is (5q > such that 

(5.3) S(k;A,B,a) e , e ,= [5„(k; A, 5)] e , e ,e lk ^ 

holds for all k G C vv/f/z Re k > arcc? Im k > (5q. 

The intuitive idea behind the proof of Theorem 15.21 is the observation that the series 
(I4.18t and (15.11) agree. However, Theorem 14 . 1 21 and Lemma l5~T1 establish convergence of 
these series in two disjoint sets of the complex plane. Therefore, to prove that both series 
define the same analytic function we perform a two-step analytic continuation invoking 
an auxiliary analytic function of two complex variables. 
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Proof. Consider the \£\ x \£\ matrix-valued function F(ki, k 2 ) with matrix elements 
F(kt, k 2 ) e>e , := [Sn(ki-,A,B)] eie ,e lk ^. 

By Theorem 15121 

F(k,k) e , e , =S(k;A5,a)e,e' 
for all k = Re k > and by Lemma l5~Tl 

(5.4) F(k,k) e>e , = S{k-A,B,a) e , e , 

for all k G C with Re k > and Im k > f3o, where (3q is defined in Lemma l5~Tl Observe 
that for any ki > the function F(ki,k 2 ) is holomorphic in k 2 G {k G C| Re k > 
0, Im k > 0}. Assume that Im k 2 > (3q with (3q defined as in Lemma l5~Tl Inspecting the 
estimates used in the proof of Lemma Bll we obtain that 

converges absolutely for all ki G C with Re ki > and < Im ki < Im k 2 + e, where 
e > is sufficiently small. Recalling ( 15.41) completes the proof. □ 

Remark 5.3. Assume that the series 

[Sn(k;A,B)] e:e ,e lk ^ 

absolutely converges in a ball B r (ko) centered at ko G C with Re k = and Im k > 0. 
Then arguments used in the proof of Theorem \5. 2\ show that 

S(k;A,B,a) eie ,= ]T [5„(k; A, £)] e , e ,e lk ^ 

for all k G i? r (ko). 

6. The Generating Function 

In this section we prove an explicit algebraic representation for the matrix- valued gen- 
erating function T(f3) defined in equation ( I2.3I) . This result is formulated below as Theo- 
rem |^1 

Let B be the canonical orfhonormal basis in C^ l+2|;r| ^K, = K, e @ JC^ /C^ +) such 
that any element h G B is uniquely associated with some edge j{h) G X U £. Moreover, 

j{h) G £ if h G K £ and j(h) G J if h G or he '. Set 

(d(j(h)) if helC s , 
v(h) = I d~(j(h)) if helC T , 
[d+(j(h)) if he!C+. 

Given a collection of matrices J\A = {M(v)} v& v we define the linear transformation 
M on the finite-dimensional Hilbert space K via its sesquilinear form 

(6.1) ^JW'WkMW " «M = ->(*.). 

I 0, otherwise. 



22 



V. KOSTRYKIN AND R. SCHRADER 



For an arbitrary > and every v £ V{Q) we set 

(6.2) A v (P):=^(I-M(v)), B v (/3) := -^(1 + M(v)). 
Define 

(6.3) A ((3):=Q)A V ((3), B(P):=QB v (0). 

vev vev 

Finally, we set 

D(P) = Z(if3;A(f3),B((3),a) 

(6.4) = I (X(iP;a) + Y(i(3;a)) - ±M (X(i/3;a) - Y(\f3;a)) 

= [I + MGH(i/3;a)]U(if3,a). 

Here the matrix Z(k; A, B, a) is denned in d3- lib , the matrices X(k; a) and Y (k; a) are 
defined in d3~T2l . U(k;a), G, and H(k;a) - in d4~15t and d4~T6l Writing the matrix M 
with respect to the orthogonal decomposition (13.31) as a 3 x 3 block-matrix 

/M n Mi 2 M 13 \ 

M = M 2 i M 22 M 23 , 

\M 3 i M 32 M 33 / 

we obtain 

fl Mi 3 e-* M i2 e-* 

I + MG a) = I + M 23 e^ M 22 e-^ 

\0 M 33 e-^ I + M 32 e- /3 ^ / 

Obviously, I + MG H(if3; a) is an entire matrix valued function in the complex variable 
(3. Moreover, 

lim det(I + MGfr(i/3;o)) = 1. 

Re/3— >+oo 

Thus, det is not identically vanishing and this in turn gives 

Lemma 6.1. The matrix valued function D{(3) is entire in (3 £ C arcc/ its determinant 
vanishes on a discrete set T> C C depending on a £ M^l an<i f/ze re? of matrices Ai = 
{M(v)} v( z\;. The set D has no accumulation points in C. In particular, the matrix inverse 
D(f3)~ 1 is a meromorphic function in (3 £ C w/f/i poles in T>. 

Now we turn to the main result of this article: 

Theorem 6.2. For a given non-compact graph Q = (V,2,£,d) with lengths a of the 
internal lines and a collection of matrices M = {M(v)} v< zy at the vertices of the graph 
the generating function T((3) defined by (I2.3t has an analytic extension toC\T> and can 
be expressed in terms of the matrix -D(/3) -1 M as follows 

ft 

(6.5) T((3) = (I 0) D{f3)~ x M I 

We turn to the proof of this theorem. First we assume that all matrices M(v) are self- 
adjoint. Let A v ((3) and B v {(3) be defined by dO . ThenA v ((3)B v ((3)^ is self-adjoint, 
since M(v) is. Observe that 

dimKer (A v (p), B v ((3)) = deg(v) - dim(Ker A v (/3) n Ker B v {j3)). 
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From dO> it follows that Ker A v (f3) n Ker B v (fl) = {0}. Therefore, the 2deg(v) x 
deg(u) matrix (A v ((3), B v {(3)) has maximal rank. Thus, the operator A(A V (@),B V (P)) 
for the single- vertex graph Q v (see Definition I4.9I > is self-adjoint. The associated scatter- 
ing matrix given by d3.15t obviously satisfies the relation 

(6.6) S(i(3;A v (f3),B v ((3))=M(v). 

Set B r (f3) = {k G C| |k-i/3| < r}. 

Lemma 6.3. The scattering matrix S(k; A v {(5), B v (fl)) is holomorphic for all 

k G (C+ \ [0,ioo)) UB r {p) 

with 

(3 

r - 



\\M(v)W 

Proof. Recalling Proposition !3.2l observe that S(k; A v ((3), B v {f3)) has a pole at k = ix, 
k G M + if and only if there is a % E such that 

5 - i) M(w) * = G + i) x ' 

that is, (J3 + >c){f3 — u)~ x is an eigenvalue of M(v). Therefore, 

|±5 < l*MI, 

which implies that the distance from the point i/3 to the closest pole of the scattering 
matrix S*(k; A v {(3), B v ([3)) is at least /3||Af(v)|| _1 . □ 

Via equations (13.61) and (I3.7t the matrices -A(/3), being defined by (16.31) define 

the self-adjoint Laplace operator A(A(f3), B(f3), a) with local boundary conditions (in 
the sense of Definition 13 .4t . 

Now we choose j3 so large that the series (I2.3t converges. Then, by (I6.6I) . the generat- 
ing function T efi i{f3) can represented in the form 

T e ,e'(/3)= 2 Sn(^;A((3),B((3)) e ^e-^, 



where the coefficients Sn{\(3;A{(3),B((3)) are defined by (14.111) and (14.121) with the 
boundary conditions d6.3t . 



Lemma 6.4. Assume that j3 > (3q with /3q satisfying (I2.12t . 77je?i ?/jere is p > .rac/z f/zaf 
f/ie series 



converges absolutely for all k G B p {(5). 

Proof. For an arbitrary e > choose p > so small that 

|^(k) ei , e2 | < ||M(«)||(1 + e) 

for all k G B p (f3), all ei, e2 G £„, and all v G V. As in the proof of Lemma BTTI for an 
arbitrary walk w G W e e '(n) we obtain the estimate 

\W(w;k)\ <m^ +1 (l + e)N +1 , 
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where 



m := max ||Af (u)||. 



In turn, this implies the bound 

\[Sn(k-,A(P),B(me,e'\<rn N+1 (l + e) N+1 \l\ N . 

\n\=N 

Therefore, 

\[Sn(k-,A(f3),B(f3))Ue'e^\ 



^E^ 1111 "^' 1 E \{Sn^A{l3),B{j3))]e,e 
\n\=N 

oo 

< ^- e -iVImka IIlillm iV+l (1 + £) Ar+l| X |iV_ 
Ar=0 

This series converges if 

(6.7) Imk > — !— log{m(l + e)\2\} . 

&min 



We claim that inequality (16.71) holds for all k € B p (/3) if e is chosen to be so small that 

(l + e )e&-0 < 1, 

and p > satisfies the inequality 

/9 <(/?-/3 )--Llog(l + £ ) 

^min 

Indeed, under these assumptions for any k G B p {[3) we have 

Imk>/3-p>/? + — log(l +e) > — log {m(l + e)| J|} . 

^min Q-min 

□ 

Proof of Theorem W2\ Assume the matrices M(u) to be self-adjoint. Lemma l6~4l and 
Remark 1531 imply that there is p > such that 

[Sn(k;A(P),B(P))] e ^e lk ^ = S{k; A{(3), B{(3)-a) 

holds for all k G B p ((3). Thus, the generating function T(fl) can be expressed in terms of 
the scattering matrix, 

(6.8) T(p)=S(ip;A((3),B((3);a). 

In turn, the scattering matrix can be calculated by means of Theorem l3.ll Obviously, 

Z(iP;A(p),B(p);a)=D(p) 

and 

A(j3)+PB(/3) = -M. 



Thus, (|63J follows from (I3TT41 
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Now we relax the assumption on the self-adjointness of the matrices M(v). Obviously, 
the r.h.s. of (I6.5t is a rational function with respect to the entries of the matrix M. Since 
det D{f3) does not vanish identically, we obtain the claim. □ 

Remark 6.5. Relation (I6.8t combined with the factorization formula for the scattering 
matrix on the graph [ 14] allows to determine the generating function T{(3) of walks on 
the graph Q in terms of the the generating functions associated with subgraphs ofQ. 

Remark 6.6. There is a direct way to establish (I6.5t . Indeed, observe that by (I6.4t one 

has 



(6.9) (I 0) .D(/3) _1 M [ 

.0, 



(I 0) [I + NLG H(i(3; a)]~ l M | 

,0; 



The matrix G performs the "jump " from one boundary vertex of an internal line to the 
other. A simple calculation shows that if He (3 > 0, then 



\GH(i(3;a) 



-Re /3 a n 



where a m j n = minaj. 

iei 



Therefore, the series expansion of [I + M G H(if3; a)] 1 converges absolutely for all (3 £ 
C with sufficiently large Re j3 > 0. The expression (I6.9t coincides with the series i2. 31) . 



This observation gives rise to the following generalization, where the penalty vector 
depends on the direction in which a given edge is traversed by a walk. Let a = {aj} ie i 
and b = {bi} ie j be two arbitrary penalty vectors. Set 



,0 









ika 











H(k;a,b) 

such that H(k; a, a) = H(k;a). Define now 

(6.10) T{0) = (I 0)\l + MGH(if3;a,b) 

For any nontrivial walk w = {e', i\, . . . , ijy, e} we set 




if the walk traverses the edge € X in the direction 
from the terminal to the initial vertex, 

if the walk traverses the edge £ 1 in the direction 
from the initial to the terminal vertex, 



where k G {1, . . . , |w| comb }. 

Using the arguments presented above one can easily prove the following statement. 

Theorem 6.7. For all (3 G C with Re (3 being sufficiently large the function T{(3) equals 
the generating function defined by the series (I2.3t with aj fc being replaced by Ci k . 
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7. Random Walks on Graphs 

In this section we define random walks on a non-compact graph Q endowed with the 
metric structure given by a penalty vector a. Assume that the matrices M(v) are stochas- 
tic, that is, all their entries are nonnegative and satisfy 

ki,k2 ~ 1 f° r an y e d§ e ^2 £ ^ U£ incident with the vertex v, 

fci 

where the sum is taken over all edges k% G T U £ incident with the vertex v. The external 
lines of the graph will be interpreted as initial or final states of the walk, the internal lines 
as intermediate states. 

Take an arbitrary external line e G £ and consider a sequence {X}^ =0 of random 
variables with values in the set X U £ determined by the following rule. Set Xq = e. 
Let vq = d(e). Choose randomly an element j\ of S(vq) with probability M(vo)j 1; e- 
Set X\ = j\. If ji G £ , then N = 2 and the sequence is completed. If j\ G T, then 
take v\ G d(j\), v\ / t>o- Choose randomly an element j% of with probability 

M(v\)j 2 j 1 and set X2 = j%. If j'2 G £, then iV = 3 and the sequence is completed. 
Otherwise proceed inductively. Finally, we obtain finite of infinite sequence of random 
variables. If N < 00, then {X}^ =0 is a walk in the sense of Section |2] We call this 
sequence a random walk on the graph Q from e G £ to e' = JTjv G £. 

The generating function of random walks from e G £ to e' G £ is defined by equation 
(I2.3I) . Obviously, it is monotone with respect to /3, 

T e , e >C£0 < T e , e ,((3') 

for /3 > (3'. If W e e / contains at least one nontrivial walk, then T e e /(/3) is strictly mono- 
tone with respect to (3, 

T e , e ,(/3) < T e>e ,((3') 

for /3 > /?'. 

Recall that the stochastic matrix M(v) is said to be regular if it is ergodic, i.e., if there 
is a natural number k such that the fc-th power M(v) k of the matrix M(v) has strictly 
positive matrix entries. 

Lemma 7.1. Let Q be a non-compact connected graph. Assume that each M{v) is a 
regular stochastic matrix. If in addition all diagonal elements of each M(v) are strictly 
positive, then all matrix elements ofT((3) are strictly positive for all sufficiently large 
[3>0. 

Proof. Connectedness of Q implies that all W efi i are non-empty. Given e and e' for 
T e ,e' (P) > to hold it is necessary and sufficient that there is at least one walk w G yV e ,e' 
with W(w) > 0. For the last condition to hold it is in turn sufficient that all matrices 
M{v) are ergodic and their diagonal elements are strictly positive. □ 

In the remainder of this section we will discuss several examples and introduce some 
mean values associated with random walks on the graph Q. These mean values are related 
to the generating function and its derivative evaluated at (3 = 0. However, for j3 > the 
generating function T e e i {(3) can be interpreted as a partition function (see, e.g., [ 23 1) with 
(3 being the inverse temperature. The role of the statistical ensemble is played here by the 
set W e e /(Ai) of all relevant walks from e' to e. 
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1. We leave it to the reader to verify that the mean length of a random walk from e' G 6 
to e G £ is given by 



(7.1) (|w|>=-^logT e ,e/C9) 



= -T e , e ,(/?)- 1 ^T e , e ,(/3) 

/3=0 a <° 



0=0 

The r.h.s. of il.lt can be calculated by means of Theorem 16.21 In the thermodynamic 
setting (i.e., for (3 > 0) the quantity 

-AlogT e , e ,(/?) = -T e , e ,(/3)- 1 ^T eje ,(/5) 

corresponds to the "mean length" at the temperature In the following examples we 
will consider probabilistic (/3 = 0) and thermodynamic (j3 > 0) means on equal ground. 

2. As another example we consider the following situation. We say that a walker entering 
a vertex v from the edge k and leaving through the edge j experiences a transition from 
k to j at v. Now fix a vertex v$ £ V and edges jo, k$ G <7„ satisfying the inequality 

M Mj Q ,k > °- We set 

(7.2) M(v ; X) jk = { ^^ahk * 3 7 * = fco 

v ^ u ' n \ M(v )jk otherwise 

with an arbitrary A > 0. Note that M(vq; X) is ergodic if M(vq) is, but of course not 
stochastic. Now replacing M(vq) by M(vq\ A) while leaving all other M(v) in the col- 
lection {M(v)} v& v unchanged, consider the matrix M(A) defined by i6.ll . Obviously, 
M(0) = M. Further, similar to (16.4b . we introduce the matrix D(f3; A) 

D(fi; A) = i (X(i/3; a) + a)) - ±M(A) (X(i/3; a) - F (ifl a)) , 
and define the generating function T(/3; A) in analogy with (I2.3t by 

T e , e ,(/3;A) = ^ ^(w;A)e- /3 H 

with 

l w lcomb 



W(w;X) = 11 [M(v k ;X)} 

k=0 

Obviously, Theorem l6.2l remains valid for T(f3; A) such that 



_(+)_(-) 

e k e k 

k=0 



(7.3) T(/3;A) = (I 0) Xy l M(X) (oj . 

Observe that if -D(/3) is invertible for a given /3 then D((3; A) is also invertible for the 
same (3 and all sufficiently small A > 0. We, obviously, have 

~~cL W<yW] ^ A=o = nv °j°' k °( w ) W ( w '>°) = n uo,jo,k (. w ) W ( w )i 

where n Vo j 0iko (w) > is the number of times a walker experiences a transition from k$ 
to jo at the vertex vq along a given walk w G W e . e /(.A/f ). 
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Consider the quantity 



(7.4) 



dX 



A=0 



d 



T e ^{(5)- l —T e ^{(5-\) 



A=0 



It is easy to verify that (n^ e ■ ko )(P) is the mean number of times a random walk from 
e' G £ to e € £ experiences a transition from ko to jo at the vertex vq. 

Using Theorem l6.2l the derivative in (17.41) can be calculated in a rather simple way: 



d_ 

dX 



A=0 



I 0)-^[D(P;Xr 1 M(X) 



A=0 



(I 0)(D(p)- l —D(fcX)\^D(p)- l M)io 



(I 0)D((3)- 1 -M(X) 



A=0 



Now set 



such that 



A=0 



M(v Q , j , fco) 



where the matrices X and F are defined in d3.12l i. Therefore, 



1 



A=0 2 



-10 



) (D{P)- l (X{^;a)-Y(ifra)) 



M(t7o,jo,feo)£>(/3) _1 M) | 

,0, 

+ (I O)Z>03)- 1 M(«o,j o ,Ak)) (o 



Thus, only the knowledge of the inverse D(f3) 1 is necessary to determine ^T(/3; A) 

Note that only one matrix element of M(t>o, jo, ko) is non- vanishing. 
The quantities 

\ n v ,j ,k 

jo 



A=0 
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are related to the mean values for the probability that the vertex vq is entered - during a 
walk from e! to e - via jo G S{vq) or left via ko G 5(t>o), respectively. Therefore, 

«fm= E 

= E «!,.>(/*) 

jo eS(«o ) 

= e <»^r.,*o>o?) 

fcoe5(^o) 

is the mean number of times the vertex vq is visited during random walks from e' to e. 
Similarly, 

( 7 -5) «fm= E «'jo, jo m 

j es(v ) 

is the mean number of times the vertex vo is entered and left through the same edge during 
a walk from e' to e. 

Assume now that for given e! G £ we have T e >e '(/3) > for all e G £. Set 

r.e'(/?) = E T ^'^)- 

Then, the value of the quantity 

a.6) <<h/3) = Ek:'>(/?)^H§ 

gives the mean number of visits at the vertex vq for random walks starting at e' G £ . 
Similarly, if for given e G £ we have T e e /(/3) > for all e' G £ we set 

T e .(/3) = E^,e'(/3)- 

The quantity 

(7-7) »(/?) = E«')(/5)^4§ 

is the mean number of visits of the vertex vq for walks ending at e G £. With 

t..((3) = e r e , C '(/?) = E T -'(^) = E T -^) 

consider the quantity 

«>(/?) = E«)(/?)|^| = D<')(^)^ 

Obviously, (n**)((3) is the mean number a random walk in W{Q) = ^ e ,e'^e,e'e£{Q) 
visits the vertex vq. Therefore, 
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is the mean number of vertices visited during a random walk. 



3. As a final example we consider the mean number (n|® ) any internal line i$ G X 
is traversed (in either direction) by a random walk from e' to e. For this replace aj by 
aj e M while keeping all other aj fixed and set 



a(io,fi) = {a,i(i , f-i)}iei with ai(i ,fJ,) 



Q>i n 6' 



«0 C ' 



if i / i , 
if i = io. 



Denote by T(/3; //) the resulting generating function. Then 



Id. 

- — logT ee i[P;fi) 
p d\i 



The derivative of the generating function with respect to \i can be calculated by means of 
Theorem 16. 21 thus yielding, 



l 



jtt=0 



-(I 0) (D(P)- 1 ((I-M)—X(if3;a(i ,ri) 



+ (I + M)^-Y(i(3;a(i ,fi)) 
dfi 



/i=0 



)D(P)- 1 M) . 



,0; 



Similar to the discussion of the mean number of vertices visited during a random walk 
we introduce the quantities 



«>(/?) 

«>(/3) 
«•)(/?) 



T e ,AP) 

T..(f3) ■ 



e'e£ 



Thus, (n* o e )((3) is the mean number of times the internal line iq G 1 is traversed by a 
random walk starting at e' G £ , (nf*}(/3) the mean number the internal line «o G X is 
traversed by a random walk ending at e G £. The quantity (n** ) (/?) is the mean number 
of times the internal line io G X is traversed by any random walk. 

Appendix. Random Walks on Vertices 

Here we will relate the customary notion of random walks on graphs (see, e.g., [2| or 
[33 1) to random walks considered in the present work. Recall that the customary notion of 
random walks on graphs is given by a Markov chain with vertices as states. The transition 
matrix P indexed by the vertices has a non-vanishing entry only if the corresponding 
vertices are adjacent. 

Consider a graph Q' = G'(V\T\ 0, d') with no external lines. Let P : V x V — ► 
M + U {0} be a nearest neighbor transition matrix, i.e., 

(A.l) J2P(v',v) = l for any veV 

v'ev 
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(we read from right to left) and P(v',v) > occurs only if v and v' are adjacent. 

Pick an arbitrary vertex in Q' which we denote by v^,. Let V Voo C V be the set of all 
vertices adjacent to v^, X Voo the set of the internal lines i G X incident with v^. For any 
i G T Voo let Vi G V Voo be the vertex adjacent to Voo by % G X, that is, 

either &{%) = {vo^Vi) or = (v^Voa) 

Now replace every edge i G T Voa by the external line e incident with the vertex v^. Denote 
the set of all external lines by £ and define the boundary operator 



d(j) 



d>(j), if jel', 

v, if j G £ . 



Thus, we have constructed a non-compact graph Q(V,l,£,d) with V = V \ and 
I = I'\ T Voo . Obviously, the degree of any vertex v G V being calculated for the graphs 
Q' and Q is equal. 

Let S(v) be the star graph of the vertex v G V, that is, the set of all edges j G 
1 U £ which are incident with the vertex v. Given a matrix P and v G V we define the 
deg(w) x deg(w) matrix M(w) with entries M(v)ij, i,j G as follows: 



< M(v) 



J P(v',v) with u' G d(i), v' ^ v, for ieS(v)\S, 
' |p(uoo,9(e)) for i = e G 5(u) nf. 

In particular, the matrix element M(v)ij is independent of j and by (IA.U 
(A.2) J^M(u)y = l for all j E £(v). 

A converse construction is also possible. Assume that a non-compact graph £ has 
£ ^ and any two vertices of the graph are adjacent by no more than one internal line. 
Further, assume that all matrix entries < M(v)ij are independent of j, that is, in each 
matrix M{v) all columns are equal, and the equality dA.2t holds. Consider the graph Q' 
without external lines obtained from Q by replacing each external line e by an internal 
incident with an additional vertex such that its vertex set V = V U {foo}- Now for 
any v', v G V we set 



(A.3) P(y',v) 



' M{v)ij for i G 2 : v',v G d(i), 

M(v)ij for i G £, v' = Voo,v = d(e), 

\£\~ x for v = Voo, v' G V Voo , 

otherwise. 



Then P is a nearest neighbor transfer matrix. 
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